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ABSTRACT 

We present parameterized broadband spectral models valid at frequencies between 30- 
300 MHz for six bright radio sources selected from the 3C survey, spread in Right Ascension 
from — 24 hours. For each source, data from the literature are compiled and tied to a common 
flux density scale. These data are then used to parameterize an analytic polynomial spectral 
calibration model. The optimal polynomial order in each case is determined using the ratio of 
the Bayesian evidence for the candidate models. Maximum likelihood parameter values for 
each model are presented, with associated errors, and the percentage error in each model as 
a function of frequency is derived. These spectral models are intended as an initial reference 
for science from the new generation of low frequency telescopes now coming on line, with 
particular emphasis on the Low Frequency Array (LOFAR). 

Key words: Radiation continuum:general - methods :observational - methods:statistical 
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1 INTRODUCTION 

In order to quantitatively combine and contrast data from inde- 
pendent telescopes and surveys, often at multiple frequencies, it 
is necessary to have a standard calibration scale to form compar- 
isons. This is especially important at frequencies below 300 MHz 
and above 15 GHz where the widely used Baars et al. (1977) ra- 
dio flux density scale is incomplete. For the new generation of low 
frequency telescopes such as the Low Frequency Array (LOFAR; 
van Haarlem et al. in prep) it is becoming increasingly necessary 
to provide a broadband spectral reference for initial science, so that 
both archival and future measurements can be quantitatively com- 
pared to these new data. In addition to an absolute scaling, such 
telescopes require a well-defined set of calibrators spread in right 
ascension (RA) to allow for quasi-simultaneous broadband calibra- 
tion of field observations. Here we present a set of parameterized 
models for six broadband calibrators covering frequencies from 
30 - 300 MHz and RAs from - 24hrs. We focus on the north- 
em sky, and in particular on the applicability to LOFAR. This set 
of calibrators forms a flux scale that will be the basis of a major ef- 
fort to develop an all-sky, broadband calibration catalog. The initial 
description given here will be continuously refined as new LOFAR 
data accumulate. 



2 CALIBRATION OF LOW FREQUENCY TELESCOPES 

Radio interferometers operating at low frequencies face a substan- 
tial calibration challenge. Strong ionospheric phase corruptions are 
common, especially below 100 MHz. For large-scale survey work 
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in particular, it is important that the processing of raw visibility data 
from the telescope can be automated. In order to jump-start such an 
automatic calibration and imaging process for any arbitrary field, a 
pre-existing model of the brightest sources in the field of view is 
required. Such a model must be intrinsically frequency dependent, 
since modern radio telescopes are inherently broadband in nature, 
with tremendous fractional bandwidths. For example, LOFAR rou- 
tinely observes from 30 — 240 MHz, and is capable of observing as 
low as 10 MHz. Over such a broad range, the flux scales must be 
tied to a well-understood set of reference sources with spectral en- 
ergy distributions that are well understood across the full bandpass. 
In the case of LOFAR, the production of such an all-sky broadband 
catalog is the key goal of the Multifrequency Snapshot Sky Survey 
(MSSS; Heald et al. in prep). 

The reference sources which form the basis of the broadband 
flux scale must be selected for suitability as high-quality calibration 
targets. Several factors are relevant. First, the source should dom- 
inate the visibility function. In addition to high flux density, sep- 
aration of contaminating flux from sources away from the point- 
ing centre ('off-beam') can be improved in two further ways (i) 
averaging in time and frequency to smear out the contributions of 
off-beam sources on longer baselines; and (ii) the "demixing" tech- 
nique (van der Tol et al. 2007), which has been adopted for use with 
LOFAR data. Secondly, the source should be compact compared 
to the angular resolution of the instrument, to allow simple mor- 
phological calibration models. Well-known sources such as Cyg A 
and Gas A have extremely complex morphologies, making calibra- 
tion of an array with arcsecond angular resolution difficult. Thirdly, 
these calibrators must be spread in right ascension (RA) to allow for 
quasi-simultaneous broadband calibration with field observations. 

With these considerations in mind we searched the 3G (Edge 
et al. 1959) and revised 3G (3GR; Bennet et al. 1962) catalogues 
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Figure 1. 3C48 SED. (a) Data fitted with Hnear (1°) model, (b) data fitted with second order model, (c) data fitted with third order model and (d) data fitted 
with fourth order model. ML parameters for each fit are listed in Table 3. Dashed lines indicate the upper and lower bounds of the LOFAR frequency band. 



for an initial list of bright compact sources, with the criteria that (1) 
they must be at dechnations greater than 30°, (2) they must have a 
flux density at 178 MHz greater than 20 Jy and (3) they must have 
an angular diameter less than 20 arcseconds (compact compared to 
the naturally weighted resolution of the Dutch LOFAR array). 

These criteria result in an initial sample of six sources, of 
which we exclude one based on other data from the literature show- 
ing more substantial extension than indicated in 3CR (3C69; Poo- 
ley & Henbest 1974) and we include one additional source based on 
other data from the literature indicating that the extension listed in 
3CR is an overestimate (3C286; e.g. Pearson et al. 1985). The final 
sample is listed in Table 1. Source extensions from 3CR are listed 
in Column [5] for each object. We note that high resolution obser- 
vations (e.g. Akujar & Garrington 1995) confirm that the source 
structure in 3C48 is on sub-arcsecond scales, whilst 3C147, 3C286 
and 3C295 have structure on scales < 5". 3C196 has two domi- 
nant components separated by about 6", as well as complex diffuse 
structure with a (precessing) jet morphology (Reid et al. 1995). The 
structure in 3C380 is known to be dominant on scales of ~ 16", 
making it the most extended object in this sample (Reid et al. 1995). 



3 FLUX SCALES 

The data used for spectrum fitting are listed in Table 2. In order to 
provide a common flux scaling, these data have been revised onto 
the flux scale of Roger, Bridle & Costain (1973; hereafter RBC) be- 
low 325 MHz. This scale has been chosen to avoid the suggested is- 
sues (e.g. Rees 1990a) with the secular decrease in the flux density 
of Cas A at low frequencies (< 100 MHz) inherent in the widely 
used Baars et al. (1977; B77) scale. 

At low radio frequencies most data are tied to the RBC or 
Kellerman, Pauliny-Toth & Williams (1969; KPW) scales. The cor- 
rection factors for moving between these scales at v < 325 MHz are 
listed in Table. 2. At v > 325 MHz the RBC and KPW scales are in 
agreement and consequently such data, where calibrated on the B77 
scale, are corrected using a polynomial fit to the correction factors 
listed in B77 onto the KPW scale. Data from WENSS (Rengelink 
et al. 1997) have been corrected using an average correction fac- 
tor to bring them onto the B77 scale and a further scaling to bring 
them onto the RBC scale. The 6C, 8C and MIYUN surveys are cal- 
ibrated on the RBC scale in their original form, and the Bologna 
survey (Colla et al. 1970) is calibrated on the KPW scale which 
is consistent with the RBC at 408 MHz. Data from Aslanyan et al. 
(1968) are scaled using the ratio of the stated flux densities for the 
calibrator sources (3C348 & 3C353) in the original paper to the pre- 
dicted values at 60 MHz from the spectral models for these sources 



Table 1. Calibration Source Sample. 



Source 



RA 

(J2000) 



Dec 

(J2000) 



(Jy) 



(arcsec) 



3C48 


0137 41.3 


+33 09 35 


50±11 


< 1 


3C147 


05 42 36.1 


+49 5107 


63 ±12 


<12 


3C196 


08 13 36.0 


+48 13 03 


66 ±20 


<12 


3C286 


13 3108.3 


+30 30 33 


21* 


<3* 


3C295 


14 1120.5 


+52 12 10 


74 ±15 


<12 


3C380 


18 29 31.8 


+48 44 46 


70±10 


<20 



' values from the revised 3C catalogue (3CR; Bennet 1962). 
unadjusted flux densities at 159 MHz ^ Pearson et al. (1985). 



in RBC. Data from Scott & Shakeshaft (1971) are corrected onto 
the scale of Artyukhetal. (1969) and then onto the RBC scale using 
the factors listed in Tables III & IV of RBC, this is subject to the 
caveat that the difference in flux densities from 81.5 to 86 MHz is 
assumed to be negligible compared to the uncertainty in these fac- 
tors (~ 3 per cent). Where applied, the scaling factors in each case 
are listed in Table 2. The original flux densities for the sources from 
the 3C catalogue (Edge et al. 1959) have not been included in the 
model fitting. The large size of the errors associated to these data is 
such that they have no influence on the parameter estimation. 

Additional data are available at 12.6 — 25 MHz from the UTR- 
1 telescope (Braude et al. 1970a;b), calibrated on the Gravoko 
scale. These data have not been used in the fitting, primarily be- 
cause the discrepancy between the Gravoko and RBC flux scales is 
not only frequency but also flux density dependent and there is no 
complete revision scale available. For a discussion see RBC. 



4 SPECTRAL MODEL 

A spectral model of the form 

log5 = logAo +Ai log V +A2log^ V H 

was used. The model was applied in linear frequency space, i.e. 



s[Jyl =^01110^''°'''''^"°"'"'': 



in order to retain Gaussian noise characteristics. Both determina- 
tion of the optimal order (A') of polynomial model and maximum 
likelihood parameter estimation were performed using a Markov 
Chain Monte Carlo (MCMC) implementation. We used a simulated 
annealing method, through the METRO algorithm (Hobson & Bald- 
win 2004), to employ a Bayesian inference approach, where Bayes' 
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Figure 2. Best fitting models for calibrator sources. Data from the literature is shown in black with the selected best-fit model overlaid as a solid line. The area 
enclosed by the shaded region indicates the flux densities allowed at each frequency by one sigma uncertainties on the parameters. The edges of the LOFAR 
bandpass (low- and high-bands) are indicated by vertical dashed lines. 



Table 2. References for data used in spectral fitting. Column [1] frequency; 
column [2] reference; column [3] connection factor applied to original data 
for conversion to RBC flux scale. 



normalizing the posterior over the prior volume, such that 



Freq. 


Ref 


factor 


10 MHz 


Bridle & Purton 1968 


1.20* 




Roger, Bridle & Costain 1973 


- 


22.25 MHz 


Roger, Costain & Lacey 1969 


1.15* 




Roger, Bridle & Costain 1973 


- 


38 MHz 


Kellerman, Pauliny-Toth & Williams 1969 


1.18* 




Rees 1990 (8C) 


- 


60 MHz 


Aslanyan et al. 1968 


1.04* 


81.5 MHz 


Scott &Shakeshaft 1971 


0.90* 


86 MHz 


Artyukhetal. 1969 


0.94* 


151MHz 


Baldwin etal. 1985 (6C) 


- 


178 MHz 


Kellerman, Pauliny-Toth & Williams 1969 


1.09* 


232 MHz 


Zhang et al. 1997 (MIYUN) 


- 


325 MHz 


Rengelink et al. 1997 (WENSS) 


0.90* 


408 MHz 


CoUa et al. 1970 


- 


750 MHz 


Kellerman, Pauliny-Toth & Williams 1969 


- 


960 MHz 


Kovalev et al. 1997 


0.96t 


1400 MHz 


Kellerman, Pauliny-Toth & Williams 1969 


- 



from RBC; ^fromB77; ^" see text for details. 



formula, 



Pr(e|£),H) = 



Pr(Z)|0,/f)Pr(0|/f) 
Pr(D|/f) 



is used to test an hypothesis, H, parameterized by using a set of 
data, D. Here Pr(0|D,//) = P{@) is the posterior probability dis- 
tribution of the parameters, Pr(£)|0,/f) = J^{&) is the likelihood 
and Pr(0|//) = 11(0) is the prior probability distribution, which in 
this case is simply used to restrict the volume of parameter space 
being sampled. 

The Bayesian evidence, Pr(£)|/f ) = Z, is a factor required for 



Z= /^(0)n(0)d'^0, 



where M is the dimensionality of the prior volume, here M = A' + 1 . 
For parameter estimation the evidence factor can be neglected as 
it is independent of the model parameters. Maximum likelihood 
(ML) or maximum a posteriori (MAP) parameter values can be 
obtained by sampling the normalized distribution in each case to 
determine the peak in parameter space. However, in model selec- 
tion the evidence becomes important for ranking different mod- 
els based on a common dataset. It can be seen from the previous 
equation that the evidence represents the average of the likelihood 
over the prior, and therefore favors models with high likelihood 
values throughout the parameter space and penalizes models with 
regions of very low likelihood. This is equivalent to numerically 
implementing Occam's razor, whereby larger evidence values are 
returned for simple models (i.e. fewer parameters) with compact 
parameter spaces, compared to more complex models - unless the 
more complex model provides a significantly better fit to the data. 
Selecting between models, say /fo and H\ , based on their evi- 
dence can be done using the ratio, 

Pr(//o|£>) ^ Pr{D\Ho)Pi{Ho) ^ Zp Pr(/fo) 
Pr(Wi |£») Pr(D|/f I )Pr(/f i ) Z, Pr(//i ) ' 

where Pr(//o)/Pr(/fi) is the ratio of prior probabilities. This ratio 
can be set before any conclusions have been drawn from the data; in 
many cases there is no reason to favor one particular model a priori 
and consequently this factor can be set to unity. In this circumstance 
the model selection can be based solely on the ratio of evidences. 

In this work, for each model, priors were assumed to be uni- 
form and separable and ML (MAP) parameters were determined 
initially using the METRO sampling algorithm. Once parameter val- 
ues had been determined, the evidence in each case, Z, was calcu- 
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lated over a ±3cT prior volume centered on the ML parameter val- 
ues, with <J determined for each parameter directly from the pos- 
terior distribution. The evidence calculation was repeated multiple 
times in each case in order to assess the variance of the evidence. 
Evidence ratios (also known as Bayes factors, or the odds) were 
then used to determine the optimal polynomial fit based on the Jef- 
freys scale (Jeffreys 1961), see § 4.1. In practice we take Aln Z > 1 
as our threshold for selecting the best model; this choice is justified 
in Section 5. 



4.1 Requirements for model selection 

The requirement to use a model of increased complexity (i.e. poly- 
nomial of higher order) depends upon the degree to which the evi- 
dence increases relative to the next lowest order, see Column [9] of 
Table 3. On the original Jeffreys scale (Jeffreys 1961) an increase 
of a factor of 3 (i.e. AlnZ ^ 3) is considered substantial evidence 
to prefer the higher order model and can be considered equivalent 
to a 99.7% confidence result. Revised versions of the Jeffreys scale 
(e.g. Gordon & Trotta 2007) divide the level of support into cate- 
gories where it is considered as either 'inconclusive' (AlnZ < 1), 
'weak' (1 ^ AlnZ ^ 2.5), 'moderate' (2.5 ^ AlnZ ^ 5), or 'strong' 
(AlnZ ^ 5). 



5 RESULTS 

The ML parameters and evidence values for each polynomial fit to 
the source spectra are listed in Table 3. An example of the dif- 
ferent orders is shown for 3C48 in Fig. 1. The poor fit of the 
linear and 2° polynomial model is evident by eye. This is also 
reflected in the values of the evidence for these models: an evi- 
dence ratio, and hence difference in the logarithm of the evidence, 
of lnZ3o — InZio > 100 indicates a definitive preference; a differ- 
ence of lnZ3o — lnZ2= = 3.16 is substantial evidence for preferring 
the 3° model above the 2° model. The fractional evidence ratio, 
lnZ4o —InZ^o = —3.55, between the 4° and 3° polynomial mod- 
els indicates that the 3° model is still preferred. In this case the 
goodness-of-fit is not diminished by the 4° model, but there is no 
evidence in the data to support the use of the extra parameter and 
hence the model is penalized. 

In general the results for this sample are easily interpreted, 
with Bayes factors of AlnZ > 3 clearly indicating a preferred order 
of polynomial in most cases. When comparing different polyno- 
mial order fits to the 3C295 and 3C380 datasets the Bayes factors 
are less conclusive than in other cases. A difference of AlnZ = 1 .6 
between the third and fourth order models in the case of 3C295 is 
intermediate to the 'weak support' category. Although the support 
for moving to the higher order model is weak, it is not inconclu- 
sive and so in the context of the work here we choose to prefer the 
fourth order model. In the case of 3C380, a value of AlnZ = 0.64 
is securely in the inconclusive category and so we prefer the lower 
order model in this instance. Best fitting spectral models for the six 
calibrator sources are shown in Fig. 2. 



5.1 Error budget 

Errors on individual parameters for each fit were determined di- 
rectly from the posterior distribution and are listed in Table 3. The 
uncertainty in the model due to these errors was derived analyti- 
cally using differential error propagation and the I a bound on the 
model in each case is illustrated in Fig. 2 as a blue shaded area. We 
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3C2a6 
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3C295 
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Figure 3. Percentage error in each model as a function of frequency from 
30 -240 MHz. 



illustrate the percentage eiTor of each model as a function of fre- 
quency from 30 - 240 MHz (the LOFAR band) in Fig. 3. It can be 
seen that from currently available data not all of the six calibrators 
are suitable for calibration at the low end of this frequency range 
(e.g. 3C147 & 3C295), where their models possess high percentage 
uncertainty. 

5.2 Notes on Individual Sources 

3C380 This source has data at 10 MHz in Bridle & Purton (1968) 
but the very low flux density (168 Jy) indicates that the spectrum 
turns over sharply below 20 MHz. The effect of this turn over is 
marginal above 30 , but would require significantly increased com- 
plexity in the model. Consequently these data have been excluded 
from the fit. 



6 CONCLUSIONS 

We have presented parameterized broadband spectral models for 
six bright radio sources selected from the 3C survey between 30- 
300 MHz, spread in Right Ascension from — 24 hours. For each 
source, data from the literature have been compiled and tied to a 
common flux density scale. These data have then been used to pa- 
rameterize an analytic polynomial spectral calibration model. The 
best fitting polynomial model order in each case has been deter- 
mined using the ratio of the Bayesian evidence for the candidate 
models. Maximum likelihood parameter values with associated er- 
rors have been presented. The percentage error in each model as a 
function of frequency has been derived and is illustrated in § 5.1. 
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Table 3. Column [1] lists the order of the polynomial fit; columns [2 — 6] list the fitted ML polynomial coefficients; column [7] lists the reduced ^ value; 
column [8] lists the natural logarithm of the evidence for a 3(7 prior volume. The selected best-fit model is high-lighted in each case. 



Order 



Ao 



As 



A4 






InZ 



A(lnZ) 



3C48: 

1° 
2° 

40 

3C147: 

2° 

40 
3C196: 

r 

3° 
40 

3C286: 



3C29S: 

2° 
3° 



5° 
3C380: 

2° 
3° 
40 



43.874 ±0.879 
62.821 ±1.642 
64.768 ±1.761 
63.910± 1.864 

60.5 17 ±1.474 



-0.349 ±0.0 11 
-0.284 ±0.0 15 
-0.387 ±0.039 
-0.394 ±0.045 

0.016 ±0.028 



-0.374 ±0.026 
-0.420 ±0.031 
-0.391 ±0.093 

-0.514± 0.046 



0.181 ±0.060_ 

0.185 ±0.075 



21.68 
1.16 

~0.014±0.118 0.15 

2.35 



-135.00±0.03 
-33.02 ±0.06 
-29.86 ±0, 

-33.41 ±0.47 

-35.28±0.14 




B6.738± 2.490 

66.494 ±1.915 

76.641 ±1.227 

|83. 084 ±1.862 

83.011 ±1.787 

83.776±2.214 



-0.022 ±0.030 

-0.041 ±0.046 

-0.719±0.012 
-0.699 ±0.014 
-0.676 ±0.029 
-0.677 ±0.033 



-1.012±0.167 

-0.952±0.109 



0.549 ±0.170 

0.625 ±0.245 



0.24 
-0.124±0.249 0.26 

2.80 



-0.110±0.02J 

-0.107 ±0.023 
-0.139±0.073 



-0.039±0.041 
-0.027 ± 0.045 



0.035 ±0.073 



0.50 
0.54 



-29.59 ±0.47 

-30. 15 ±0.48 

-47.25 ±0.08 
-36.89 ±0.09 
-38.16±0.14 
-40.37 ±0.50 



97.489 ±2. 177 
100.950 ±2.454 
197.763 ±2.787 



-0.347±0.016 
-0.517±0.035 
-0.582 ±0.045 



-0.362 ±0.028 
-0.497±0.041 
-0.298 ±0.085 



0.360 ±0.066 
0.583±0.116 



-0.363±0.137 



75.682± 1.537 
75.233 ±1.483 
74.386 ±1.595 



-0.767 i0.013q 

-0.772 ±0.0 12 
-0.788±0.033 
-0.787 ±0.034 



0.039 ±0.021 
0.041 ±0.020 
0.104 ±0.067 



0.024 ±0.047 
0.030 ±0.051 



1.02 

1.11 

-0.084±0.082 1.14 



-3«. Hi 0.06 

-37.47 ±0.10 
-39.22 ±0.15 
-40.96 ±0.40 



5.69 

-0.59 



10.36 

-1.27 
-2.21 



1° 

2° 


27.893 ±0.653 
28.230 ±0.708 


-0.258±0.017 
-0.208 ±0.035 


-0.077 ±0.045 


— 


_ 


1.46 
1.32 


-33.93 ±0.08 
-34.09 ±0.14 


-0.16 




^7.477 ±0.746 


-0.158±0.033 


0.032 ±0.043 


-0.180±0.052 


^^^^^^^1 




t30.51±0.21 


3.58 


40 


27.591 ±0.911 


-0.144 ±0.038 


0.005 ±0.097 


-0.187±0.054 


0.021 ±0.086 


0.48 


-32.58 ±0.49 


-2.07 





These spectral models are intended as an initial reference for sci- 
ence quality data from the new generation of low frequency tele- 
scopes, such as LOFAR, now coming on line. In this context we 
have shown that two of these sources lead to unacceptably high flux 
scale uncertainty at frequencies below 70 MHz (3C147 Sl 3C295), 
and we also note that 3C380 may be unsuitable for precision cali- 
bration at higher frequencies where its angular extent becomes an 
issue. 
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